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We investigate the spin-exciton modes in the superconducting d x 2_ y 2 state of CeMIn 5 heavy 
fermion compounds found at the antiferromagnetic wave vector by inelastic neutron scattering. We 
present a theoretical model that explains the field dependence for both field directions. We show 
that the recently observed splitting of the spin exciton doublet in CeCoIns into two non-degenerate 
modes for in-plane field appears naturally in this model. This is due to the spin anisotropy of g- 
factors and quasiparticle interactions which lead to different resonant conditions for the dynamic 
susceptibility components. We predict that the splitting of the spin resonance doublet becomes 
strongly nonlinear for larger fields when the energy of both split components decreases. For field 
along the tetragonal axis no splitting but only a broadening of the resonance is found in agreement 
with experiment. 

PACS numbers: 74.25. Jb, 71.27.+a,72.15.Qm 



I. INTRODUCTION 

In unconventional superconductors quasiparticle exci- 
tations that determine low temperature thermodynamics 
and response exhibit a generally anisotropic gap A(k) 
with posible node lines on the Fermi surface (FS). In 
addition to these single particle excitations collective 
excitations may appear. The collective oscillation of the 
superfluid density or 'Higgs mode' that belongs to the 
singlet spin sector is diffcult to observe. In addition 
collective spin triplet excitations may be present below 
the gap edge that are formed as bound states due to 
quasiparticle interactions. They are accessible directly 
by inelastic neutron scattering (INS) and have been 
found in a considerable number in unconventional 
superconductors. In heavy fermion compounds their 
typical energy is in the range of just one mcV. The 
most clearcut example belongs to the class of 115 
superconductors, CeMIns (M= Rh, Ir and Co). They 
have attracted great interest because of coexisting and 
competing antiferromagnetic (AF) and superconducting 
(SC) states^"— and in particular due to the possible 
existence of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
phase in CeCoIn^. 

The highest superconducting critical temperature of 
this family appears in CeCoIii5 with T c = 2.3K, and 
the question of the gap symmetry in these compounds 
has been intensely discussecP^. A powerful indirect 
method to study the unconventional gap symmetry is 
provided by the spin resonance peak which may appear 
in inelastic neutron scattering (INS). Such a pronounced 
spin resonance has been observed in CeCoIn 5 and La 
substituted crystals at w r /2Ai = 0.65 in the supercon- 
ducting state by INS where 2Ai is the main quasipar- 
ticle gap obtained from tunneling experiments^. It is 
confined to a narrow region around the AF wave vec- 
tor Q = (5,5, ^Jp^- A theoretical calculation of the 
dynamical spin response using realistic Fermi surface^ 



shows that the resonance can appear only for the d x i_ y i- 
wave gap symmetry but not for d xy - type gap because the 
precondition A(k + Q) = — A(k) for a spin resonance^ 
is only fulfilled in the former case. Indeed the d x 2_ y 2 gap 
symmetry has also been found by thermal conductivity 7 
and specific heat measurements^ in rotating magnetic 
fields. 

The existence of a spin exciton resonance in the SC 
phase is a well established many-body effect which has 
also been observed in other unconventional supercon- 
ductors like high-Tc cuprat es^, heavy fermion metals 
UPd 2 AlJ2E£l and CeCuaSiJ™, and in particular in 
many Fe-pnictide compounds^. The appearance of the 
resonance depends sensitively on the type of unconven- 
tional Cooper-pairing and provides a powerful criterion 
to eliminate certain forms of pairing when a resonance 
is observed. These interpretations, however, all depend 
considerably on theoretical phenomenological model fea- 
tures like FS nesting properties, nodal positions and mo- 
mentum dependence of the gap as well as the size and 
anisotropy of quasiparticle interactions. 

The spin exciton is a triplet excitation in the isotropic 
case since it should appear as a pole or resonance in all 
three components of the susceptibility tensor. In princi- 
ple, in presence of the magnetic field it should split into 
three modes with different polarization (left and right 
handed as well as longitudinal). The field- induced split- 
ting of spin excitons has been predicted for the cuprates^ 
but was sofar never identified experimentally in any un- 
conventional superconductor. 

Recently an apparent field induced spin-exciton split- 
ting has been found for the first time in CeCoIns by 
INS and has given rise to an interesting debate. In 
Ref. [21] a splitting into two (rather than the expected 
three) modes was found for field in the tetragonal plane 
and no splitting for perpendicular field. On the other 
hand no splitting was reported in Refs. (|22l and l23|) for 
both field directions and only an increased broadening 
of the peak is observed^. Furthermore it has been pro- 
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FIG. 1. (Color online) Fermi Surface (a) and quasiparticle 
density of states (b) for zero magnetic field, here 2Ai = 
0.56Ao and 2A2 = 1.5Ao. (c) and (d) show Fermi surfaces 
for spin up and spin down, respectively, in presence of the 
magnetic field Bo z, corresponding to Zeeman energy splitting 
hg Q — g^x^BBo — 0.3Aq which is also used in subsequent 
figures. 



posed that the FFLO-type inhomogeneous superconduct- 
ing 'Q-phase'SSllSZlis due to the condensation of the lower 
spli t-off b ranch of spin excitons close to the upper critical 
field™. 

In this paper we present a theoretical analysis to ad- 
dress the field splitting of the spin resonance which have 
been observed recently in INS experiments^ for the first 
time in an unconventional superconductor. We investi- 
gate whether the conventional picture of this mode as 
a triplet bound state of quasiparticles can explain these 
observations. We clarify the appropriate conditions for 
the splitting to occur for the different field directions as 
well as the number and field dependence of the split spin 
cxciton modes. 



II. THEORETICAL MODEL 

The Anderson lattice model provides a convenient 
frame to model the heavy quasiparticle bands in CeMIn 5 
realisticalljEH. The electronic structure of CeMIns has 
been investigated by using tight binding models with ef- 
fective hybrid isation for f-electrons and c- (conduction) 
electronsPHni. The crystalline electric field (CEF)- 
splitting is approximately three times larger than the 



quasiparticle band width (W « 4 meV), then one may 
restrict the model to consider the lowest Ty 1 ' Kramers 
doublet of 4f states^. This can be described by a pseudo- 
spin a degree of freedom. Choosing its quantization 
axis along the z direction which is defined by the mag- 
netic field, the Anderson lattice model Hamiltonian for 
the two hybridized conduction and localized orbitals (c,f) 
which are doubly Kramers degenerate is given by 

« = E e k<T c kcr c ko' + eUlho + V k (cL/W + h.c.) 

ka 

+E f/ ///k t /kt/^M- (i) 

kk' 

Where c ka creates an electron with spin a in the conduc- 
tion orbital with wave vector k = (k x ,k y ,k z ). Further- 
more, e Wa = — T~L% an d sjLr = e k — > wnere e k an d e k 
are effective tight binding dispersions of the conduction 
band and the renormalized dispersion for the / band, re- 
spectively. The Zeeman splittings of bands due to the 
effective molecular fields are given by W B — h c B a z with 
h B = \9 c X^bB and W B = h B a z with h B = \glx^sB. 
Here a z is the Pauli matrix, B is the magnetic field, g c is 
the g-factor for c-electrons, g^ a the g-factor for /-electrons 
in direction a =|| , _L with respect to the tetragonal plane, 
X the Stoner enhancement factor of the homogeneous sus- 
ceptibility due to quasiparticle interactions and \xb is the 
Bohr magneton. The anisotropy of the g-factors is ob- 
tained from that of magnetization or spin susceptibilitjEH 
as gj_ /gl = 2.3 assuming that the f-electron contribu- 
tion in the magnetization dominates. For that reason we 
chose a small g c /g^ = 0.2. Furthermore creates the 
f-electron with momentum k and pseudo spin er, and 
Uff is its on-site Coulomb repulsion. Finally Vk is the 
hybridization energy between the lowest 4f doublet and 
conduction bands which contains implicitly the effect of 
spin orbit and the CEF term and is taken as momentum 
independent. 

It is knowrP3 that in the limit of Uff — > 00 where dou- 
ble occupation of the f-states are excluded, an au xiliar y 
boson defines the mean field (MF) Hamiltonian a d 33 * 34 ! 

« = E E iA^° + C/Liw + v k (4 CT / kCT + h.c.) 

ka 

+A(r 2 -1). (2) 

The MF Hamiltonian can be diagonalized using the uni- 
tary transformation, 

fka = M+,kcra+,kcr + U- t koOt-,ka 

Cka = M-.ktra + .kCT - W+,ko-a-,kCT (3) 

and as a result one can find the MF quasiparticle Hamil- 
tonian as 

h mf = £ E^ai M a ±M , +X(r 2 1), (4) 

±,kcr 
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FIG. 2. (Color online) Individual compounds (xo (Q,u/)) and 
total (xo° t (Qi tJ )) of the unperturbed pseudo spin susceptibil- 
ity in absence (a) and at presence of the magnetic field along 
the z direction (b), and for magnetic field along the x direc- 
tion (c). (Bo is the same as Fig[T|. 



where the two pairs of quasiparticle bands (pairwise de- 
generate for zero field) are given by 



E t = - 2 



'kcr 



, (5) 



and the quasiparticle mixing amplitudes are obtained 



0.8 
0.6 
Q 0.4 
C 0.2 

O 

x 0. 
-0.2 
-0.4 







■ lm X m , B = 
□ Re;f' 01 , B = 
• lmx m , B = Bo t 




O Re*", B = B i 
A \mx'°\ B = B i I 
A ReA- 1 "', B = B St 











w(A ) 



FIG. 3. (Color online) The total unperturbed physical sus- 
ceptibility (including matrix elements mn,m±; x tot (Q>u) = 
X zz + |[x +_ + X _+ D f° r zero field and finite field along z, and 
x directions. (Bo is the same as FigQ. 
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Using the parameters defined in Ref. for the above 
quasiparticle band structure wc plot the corresponding 
FS in Figflja) in the absence of a magnetic field. Since 
the Fermi level is located in the lower band, E^ a , we 
can neglect the upper band, for discussing the low 
energy spin excitations. 

The superconducting pairs are then formed from the 
heavy quasiparticles of the lower band leading to a pair- 
ing potential in the Hamiltonian according to 



u 



sc 



h.c. 



(7) 



where A k is the superconducting d-wave gap function for 
CeMIns given by 



A 



(cos k x — cos k y ). 



(8) 



By defining the new Nambu spinors as = 
the effective Hamiltonian can be written as 

/eS ^ 



here 
and 
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FIG. 4. (Color online) RPA susceptibility in absence and 
presence of the magnetic field: along the z direction (a), and 
along the x direction (b). Parameters corresponding to the 
behaviour as observed in CeCoInswith gf /g^ = 2.3 and Jq = 
2.6A ; Jq = 13.5A . (B is the same as Fig[l]). 



The propagator matrix of the conduction electrons in 
terms of the Nambu spinor in the Matsubara representa- 
tion is obtained as G(k, r) = — (Tipk(T)ip^(O)) . Using a 
standard equation of motion method one can find that 



G(k,Lu n ) = (iuj n - /3 k ) 
Explicitly written, 
G(k,r) 



Gf(k,r) Gf{k,r) 



(11) 



(12) 



where 6f (k,r) = -<T> ?k (r)^, k (0)) and = 1,2 
denote the spinor components for ±k below Eq. Q9p . 



A. Magnetic susceptibility 

The noninteracting or bare dynamical f-electron sus- 
ceptibility is defined by 



xS'(r) = -f(r)<Ti4(r)£ q (0)) > 



(13) 



0.7 



< 

3 



0.6 



0.5 





(a) 


CeCoIns ^^'^'^ 






B\\z 





0.2 



0.4 

h f B (A ) 



0.6 



0.7 



3 



0.6 



0.5 



(b) 



B\\x 



CeCoIn 5 



0.2 



0.4 

h f B (A ) 



0.6 



FIG. 5. (Color online) RPA susceptibility peak positions 
corresponding to spin resonance energies as function of field 
strength: along the z direction (a), and along the x direc- 
tion (b). Anisotropic g- factors with gf /g\ — 2.3 and quasi- 

particle interaction parameters Jq = 2.6Ao; Jq = 13.5Ao 
leading to similar behaviour as in CeCoIns for both field di- 
rections. Here the maximum h B = 0.7Ao corresponds to a 
field B = 39.7 /(g l f x) in Tesla for direction / =J_, ||. 



where 



Jq 



kcr o 



r M cra'ha' 



The hybridizing 4f states are in a 11^ ) Kramers doublet 
state and therefore their physical moment operator (in 
units of /ib) can be presented by a pseudo spin matrix a 
according to (I = +,—, z) 



(14) 



Mi = miai 



(15) 



Here mi are matrix elements and their anisotropy (m± = 
mn;m z = mi) may be directly obtained from that of 
the experimental spin susceptibility of quasiparticles^ 
according to m±/m^ = .g^/.gjj = (xl/xjj)^ = 2-3. Their 
absolute value is not significant as it only sets the overall 
scale. Then, the bare physical moment susceptibility can 
be expressed as 



x"'(q,w) = m ; m/'Xo'(q,w), 



(16) 
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where Xo (q> w ) is the pseudo spin susceptibility defined 
by 

Xo (<1: W ) = ^ Cr L' Cr l l( Ti W -,k+q ( TW-.k+q<T;W-,k ( T 1 U-,k ( T' 
M<t} 



dio 



(17) 



With the definition of the new basis as ip^ = 
(&+,ik>&+,2k>&-,2k)6-,ik) an d applying the Bogoliubov 
transformation given by 

«-,kt = u +,lk & +,lk + W-,lk&-,lk 
"-,-kJ. = -«-,lk&+,lk + «+,lk&-,lk 
a -,kl = w +,2k^+,2k - «-,2k&-,2k 
"-,-kt = V-,2k6+,2k + V+,2k&- : 2k (18) 

the effective Hamiltonian in Eq. ([9| is diagonalized by 
= PkAcPk, where Pk is a 4 x 4 matrix, composed 
of the eigenvectors of /3k- Here /?k is the diagonal matrix 
constructed from the corresponding eigenvalues, 

E tk = E f_ k = 2 



E^~EZ n ±^{E^ + EZ^) 2 + 4A k 

(19) 



and P k 1 is the matrix inverse of Pk- Furthermore, 



(^k t + ^) 2 + 4A, 2 



(20) 



Then G d (k, uj n ) = (iu n — is diagonal and we can 

write 

4 

Gll^UJn) = ]T 7 ^ lS <<%(k,U, n ), (21) 
s' = l 

where r y^ aiS , = P^ s ,Pks'a 1 - Thus, the pseudo-spin sus- 
ceptibility can now be expressed as 

k,{a} { S } 

/(^ +q ,, 2 )-/(^, 4 ) 



u -,k+qcr U- .k+qo-J u - ,kcr 1 M - ,ka' 



k+q,s 2 



(22) 



where /(e) is the Fermi function. The combination of 
prefactors in the sum of Eq. ( 22 1 are the combined co- 



herence factors arising from the hybridisation and the 
superconducting state. Finally the cartesian dynamic 
magnetic susceptibility tensor in random phase approxi- 
mation (RPA) has the form 



where J(q) is the effective quasiparticle interaction ma- 
trix with non-zero elements: J^ z = J^, j£~ = J~ + = 

jjj. Similar to the g- factors they may be anisotropic. 
Furthermore their momentum dependence can be mod- 
eled by a Lorentzian function which is peaked at the 
wave vector Q associated with the indirect hybridization 
gap^. Finally the total dynamic magnetic susceptibility 
is obtained as the trace of the tensor according to 

Xrpa{<1,u) = Xrpa(<1>u)+2 [x~rpa(<1>u) + Xrpa^^)] 

(24) 

For momentum transfer q along Q = (|, |, |) direction 
this is directly proportional to the dynamical structure 
factor S(q, lo) observed in INS experiments. 



III. NUMERICAL RESULTS 

In this section, we evaluate numerically the dynamic 
magnetic susceptibility for our model of 115 systems. We 
consider two different cases, namely magnetic field along 
z- direction (out-of-plane: c-axis) and along x- direction 
(in-plane: a-axis). For clarity we give first the relevant 
physical parameters of CeCoIns to simplify comparison 
with experimental results. For the main (lower) tunnel- 
ing gap (Fig. [TJ we have 2Ai = 0.56A = 0.92me"V® r 
A = 1.64meV for the gap amplitude in Eq. ([8]). The 
spin resonance energy for zero field is io r — 0.6 meV^or 
w r /2Ai ~ 0.65. 



ponents Xo (Qi w ) °f the pseudo-spin susceptibility ten- 
sor and its trace Xo°*(Q,w) = x§ z (Q, + \ [xt (Q, w) + 
Xo +(Q,w)]. For B -> Figga shows that \ X t~ > \xt " 
and X q z become identical, independent of the field di- 
rection due to isotropy in pseudo-spin space. Therefore 
X ot for B — > in Fig{2]a is simply three times each 
of these components. This means that the anisotropy 
enters only through the matrix elements in the physi- 
cal moment susceptibility of Eq. ( 16 1 discussed below. 



First, using Eq.(22), we calculate the individual com- 



Xrpa(<1, w) = [1-J (q)x(q, w)] x x(q, w), 



(23) 



The pseudo-spin susceptibility in presence of a magnetic 
field along the z-direction is shown in Fig(2jb and for x- 
direction in Fig(2jc. We note here that for calculating 
the pseudo spin susceptibility with magnetic field along 
the i-direction, we have applied a 7r/2 rotation along the 
y-direction (x — ¥ —z and z — > x) and similar in k-space. 
For finite field and both directions the individual com- 
ponents for each field direction start to differ because of 
the polarization of quasiparticle bands in the field. This 
leads to a change of the nesting conditions of the FS in 
presence of a magnetic field, shown in Fig.([T]b), which is 
different for each susceptibility component. Most impor- 
tantly the energy dependence of Xo~~ and Xo + shows a 
splitting in opposite manner. 

We also calculate the total bare physical suscep- 
tibility, x* 0< (Qj w ) resulting from Eq.([l6| in Fig{3] 
The anisotropy enters only through the matrix ele- 
ments m\\,m±_ in the physical moment susceptibility of 
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Eq. (16 1. It is shown at zero magnetic field (B 0), 



and also for magnetic field along z- and i-directions 
(in original axis notation), respectively. We note that 
for B — > the ^x" 1 an d \x h components are still 
equivalent as in Figj2ja but differ from \ zz due to the 
different magnetic moment for in plane and out-of-plane 
direction (see Eqfl5]). 



Finally, from Eq.([24j) the interacting RPA susceptibil- 
ity and the spectrum of excitations given by its imaginary 
part can be calculated. In general form, we obtain 



XflPA(q,w) 



m 2 zX o z 



-Xo + ) 



{X x m 2 y + X y m 2 x )xo X ~ 



4-(\ x + \ v )(xJ +Xa + ) + ^XyXa X ~ 



where the interaction parameters are defined as 

A, 



m^Jq (here / 



= x — y —\\,z =_L for the orig- 
z = y =||, a; =_L for the rotated 



inal axes, and I 
axes). If the resonance condition should be satisfied for 



both transverse parts in Eq. ( 24 1 then one must have 
myJq 1 Rj m^Jq. Since the CEF states and hence the 
anisotropy may change as an effect of substitutions in 
the pure 115 compounds the resonance signatures may 
also change accordingly. As a result of the sign change 
of the superconducting gap function (Ak+Q = —At) 
at the antiferromagnetic momentum Q the spectral 
function lmxo{Q,u) remains zero for the low frequencies 
and then shows a discontinuous jump at the onset 
frequency of the particle-hole (p-h) continuum, i.e., close 
to uj c = min(Aic+Q + Ak). This is around Ao ~ 2Ai 
, where 2A X is the gap in the superconducting DOS in 
(Figfljb) which is observed in the tunneling spectrum 
of CeCoIn 5 9 . The resonance may appear for energies 
u> < u} c , under the condition that (i) J q Rex 0q (u;) = 1 



and (ii) Imx" q ( 



(«' = +-,— 



We begin our discussion of numerical results by consid- 
ering the strongly anisotropic case with Jq <C Jq. Only 
in this case, the resonance condition is satisfied for both 
Xqq and Xoq components and not for Xoq: i- e -j a reso- 
nance doublet is possible as observed in experiment. At 
zero magnetic field a single sharp peak for the degener- 
ate resonance is observed which is shown in Fig|4] By 
applying the magnetic field the bare Xo~'Xo + suscep- 
tibilities start to split into an upper and lower branch 
(Figs. [2jb,c), and as a result the two resonance peaks 
of the doublet are revealed in the RPA frequency spec- 
trum. Because of the anisotropic interaction, these peaks 
are completely distinguishable for in-plane magnetic field 
showing a linear Zeeman splitting for small fields. On the 
other hand they merge together for the out-of-plane field 
and a single broadened peak appears whose width in- 
creases with field. Using the procedure described above 
for various fields we obtain the peak positions of the RPA 



spectrum versus magnetic field strength in FigjEJ For 
field oriented along the z-direction (Fig{5ja) we always 
have a single peak with larger broadening, and by in- 
creasing the magnetic field the peak moves to lower en- 
ergies. But for the magnetic field applied in i-direction 
(Fig{5]b) the RPA result shows always two peaks with 
narrower line widths. When the field is increased these 
peaks show first a linear Zeeman splitting and one of 
them moves to larger and the other one to lower ener- 
gies. Finally at large field nonlinear behaviour sets in 
and both start to move to lower energies. These results 
are in complete in qualitative agreement with the exper- 
imental observation for CeCoInJ^H at lower fields which 
have sofar been used only. It is clear that the linear split- 
ting observed must be modified when larger fields closer 
to the upper critical fields are applied. The absolute field 

39.7/( 5 < x) m 



scale at the maximum hL = 0.7A is B 



(25) 

Tesla for direction I =||, JL. From a comparison with the 



experimental linear splitting region up to 6T <C H^ 2 and 

the one in Fig 5a with h? B 0.35 one gets the param- 
I — 

eter g'^x = 3.3. We note that the reference scale Ao is 
only constant for fields small to the upper critical field 
{h\ 2 = 11.9 T, Hjr 2 = 4.95 T). For larger fields one has 
to scale with Ao(B) which vanishes at the upper criti- 
cal fields H^ L where the resonance energies ui r (B) also 
have to vanish for both field directions. When the curves 
in Fig.[5^,b are multiplied by the scaling function A (B) 
the field dependence ui r (B) is obtained in absolute (meV) 
units. 

Finally we mention here that our results within the 
phenomenological RPA theory depend on the model pa- 
rameters, i.e. anisotropic matrix elements and quasi- 
particle interaction energies. Changing these parameters 
leads to other interesting regimes, where the anisotropies 
of magnetic moment and quasiparticle interaction along 
the tetragonal axes play the decisive role. The RPA treat- 
ment shows that even in absence of the magnetic field an 
additional peak from the Xo q component may appear for 
suitable parameters and therefore in principle two peaks 
may exist in the RPA spectrum even at zero field. In the 
presence of the magnetic field like previous cases the de- 
generate (doublet) peak of the Xoq" an d Xo~q^ components 
splits and finally leading to three different peaks at finite 
field. We should stress here that by moving to lower in- 
teraction energy the resonance occurs just for Xq^ an d 
the splitting of the resonance vanishes. This could give 
an insight into the challenging alternative experiments 
where the splitting behavior for in-plane field has not 
been seerP^l 



IV. CONCLUSION 

We have given an explanation of the observed field 
splitting of feedback spin resonance excitations recently 
observed for the first time in the unconventional super- 
conductor CeCoIns. Our calculations are based on a RPA 
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model with tetragonal anisotropy for collective spin ex- 
citations that may also be relevant for other members of 
the 115 family. Using the appropriate hybridized bands 
and associated Fermi surface as well as the proper d x 2_ y 2 
superconducting gap symmetry the spin resonance ap- 
pears at u/2Ai = 0.73, close to the experimental value. 
Here 2Ax is the main tunneling gap indicated in the DOS 
of Fig[T]b and found in Ref. GO 

In the isotropic case the resonance is a spin triplet ex- 
citation that should split into three modes in a field as 
predicted in the case of a cuprate modeP4 However, be- 
cause of the presence of strong CEF and hybridization 
induced anisotropics of g-factors and interactions of 4f- 
type quasiparticles in 115 compounds the resonance con- 
dition may in this case not be fulfilled for all three triplet 
components. For sufficiently strong anisotropy the zero 
field resonance is only twofold degenerate transverse dou- 
blet state because the longitudinal component does not 
satisfy the resonance condition . For field in the tetrag- 



onal plane the doublet splits with a linear Zeeman effect 
for low fields, one branch moving to higher the other two 
lower energies. However, for larger fields a crossover to 
nonlinear field dependence with both split resonance en- 
ergies decreasing sets in. For field perpendicular to the 
tetragonal plane no splitting but only a broadening of 
the doublet resonance excitation appears. These salient 
features of our model calculation correspond closely to 
the experimental observations in CeCoIns performed for 
small fields^. It would be very interesting to investigate 
the predicted nonlinear field dependence for larger fields. 

Finally we note that other scenarios are possible de- 
pending on the anisotropics and strengths of quasipar- 
ticle interactions where for example the in-plane split- 
ting of the resonance disappears and a single peak with 
anomalous broadening as for the out-of-plane field is ob- 
served. In fact this behaviour was proposed in an al- 
ternative experiment!^ and further experimental as well 
as theoretical investigations are necessary to clarify this 
issue. 
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